Effects of surface forces and phonon dissipation in a three-terminal nanorelay 
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We have performed a theoretical analysis of the operational characteristics of a carbon-nanotube- 
based three-terminal nanorelay. We show that short range and van der Waals forces have a significant 
impact on the characteristics of the relay and introduce design constraints. We also investigate the 
effects of dissipation due to phonon excitation in the drain contact, which changes the switching 
time scales of the system, decreasing the longest time scale by two orders of magnitude. We show 
that the nanorelay can be used as a memory element and investigate the dynamics and properties 
of such a device. 
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I. INTRODUCTION 



Nanoelectromechanical systems (NEMS) is a rapidly 
expanding research field with high potential for 
applications. 1 Typically, the length scale of these systems 
is in the nanometer range and there is a strong coupling 
between the mechanical and electrical degrees of free- 
dom. Key properties of NEMS are high natural frequen- 
cies, operation at ultra-low power, and low dissipation. 
Carbon nanotubes (CNTs) have extraordinary properties 
that make them ideal building blocks of NEMS. They are 
light, have very high Young's moduli^ and can be elas- 
tically deformed without breaking^ In addition, CNTs 
have a wide range of electric properties: they can for in- 
stance be either semi-conducting or metallic depending 
on detailed structured 

In this article, we analyze a nanorelay system in which 
a conducting multi-wall CNT (MWNT) is placed on a 
terrace in a silicon dioxide substrate and connected to 
three electrodes. This system has previously been studied 
by Kinaret et. al^ who introduced a simple model and 
discussed the characteristics of the device. In the present 
work, we focus on the effects of van der Waals and short 
range forces in the system. We show that the magnitudes 
of these forces are comparable to those of the elastic and 
electrostatic forces and that they need to be included in 
a more realistic model of the system^ A brief account 
on this topic has been presented in Ref. 0. The surface 
forces alter the system's behavior, and the characteristics 
prove to be useful for an application of the relay as a 
memory element. 

We also investigate the dynamic performance of the 
proposed memory element and estimate its write times. 
We show that the switching time depend sensitively on 
dissipative surface processes related to tube-surface in- 
teractions when the tube impinges on the drain contact, 
and analyze this in more detail using a simple model 
phonon emission. The dissipation associated with this 
interaction changes the switching time scale of the relay, 
and consequently, change the typical write time of the 
proposed memory element. 



II. THE MODEL 



A schematic picture of the system is presented in 
Fig. 2] A metallic multi-wall carbon nanotube of length 
L is placed on a terraced silicon dioxide substrate such 
that its free end is distance h from the drain contact when 
the tube is horizontal. The fixed end of the nanotube is 
connected to a source electrode, to which the electric con- 
tact is purely ohmic with impedance Z. Two additional 
electrodes (gate and drain) are positioned on the lower 
level. The position along the horizontal axis is labeled z. 
The drain electrode is positioned at z^ — L and the gate 
at z g < L. The ratio z g /zd — z r gives the relative posi- 
tion of the gate. The free end of the tube can move and 
the electrodes are used to control its motion. The deflec- 
tion of the nanotube tip from the horizontal position is 
labeled x. Typical length scales of the system are tube 
length L = 50 — 150 nm, and step height h = 5 — 10 nm. 
By applying a voltage to the gate and/or the drain elec- 
trodes, an excess charge q is induced in the tube. This 
causes an electrostatic force that deflects the tip of the 
tube towards the drain electrode. 

The tube is modeled using classical continuum elastic- 
ity theory^ considering the CNT to be an elastic can- 
tilever with only the lowest vibrational eigenmode ex- 
cited. Furthermore, we assume that the bending profile 
of the tube affected by an external force F is the same 
as for free oscillations. Then the potential energy of the 
tube can be expressed in terms of tip deflection x accord- 
ing to V — kx 2 /2 where the effective spring constant is 
given by k « 3EI/L 3 , where I = n(D^ - £>f)/64 is the 
moment of inertia of the cantilever, and E is the Young's 
modulus, experimentally measured to be of the order of 
1 TPa.SiiS, The equation of motion for the tip of the tube 
is that of a forced harmonic oscillator with an effective 
mass given by m c s = fc/fi 2 « 3Mcnt/(1-875) 2 where f2 
is the lowest vibrational eigenfrequency and Mqnt is the 
total mass of the CNT-cantilever. In addition, a phe- 
nomelogical viscous damping j^x is introduced into the 
equation of motion, 7d is estimated from Q-factors from 
experimental data. 11 Finally we arrive at the equation of 
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Figure 1: (a) A schematic picture of the system. One end of a conducting multi-wall carbon nanotube of length L, inner 
diameter Di, outer diameter D a is attached to a source (s) electrode above a step of height ft in a silicon dioxide substrate. 
The other end is free to move. The displacement x of the nanotube tip is measured towards the substrate. A gate (g) and a 
drain (d) electrode are placed beneath the tube and are used to control the motion of the tip. (b) An equivalent circuit for 
the system in Fig. 1(a). The impedance Z is ohmic and describes the tube-source coupling. The tube-gate coupling is purely 
capacitive with the capacitance C g (x) whereas the tube-drain coupling is a tunnel junction with tunneling resistance Rt(x) 
and capacitance C'd(x). V s , g ,d are the electrostatic potentials on the electrodes and we choose Vd = 0. 
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m cS x 



-kx - 7 d ± + F c + F vdw + F m 



(1) 



where F c is the capacitive force, -FVdw is the van der 
Waals force, and F sr are the short range forces due to 
wave function overlap. These different forces will be dis- 
cussed in more detail in the next section. The charge 
dynamics in the tube is governed by 
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where V s and V g are the electrostatic potentials on the 
source and the gate, respectively, Z is the source junction 
impedance, C g and C d are the capacitances between the 
tube and the gate and the drain respectively, and / s d 
is a stochastic tunneling current between the tube and 
the drain electrode. Due to the much larger tube-gate 
separation tunneling to the gate is neglected. We choose 
Vd = Q. 

The tunneling rates of the stochastic tunneling current 
are calculated using the theory of Coulomb blockade in 
the presence of an electromagnetic environment. Follow- 
ing Ingold and Nazarovj^S, assuming an ohmic environ- 
ment and zero temperature, the tunneling rate T for a 
given deflection x is determined by 
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where R? is the tunneling resistance and V is the electro- 
static potential difference between the tube and the drain 
contact (V ~ V s ). Furthermore, P{E) is the probability 



for energy exchange between the tunneling electron and 
the environment and is self-consistently determined from 
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where E c (x) — e 2 /[2Cs(.x)], where C^(x) = Cd + C g is 
the total capacitance between the tube and the contacts 
and i?K = 2irh/e 2 « 25.8 kO is the von Klitzing constant. 
The tunneling resistance is approximated by 

R T = Roe<- h - x V x , (5) 

where Rq is the tunneling resistance when the tube is in 
mechanical contact with the electrode, A is the material 
dependent tunneling length of the order of 0.5 A and 
h is the tip-contact separation at zero deflection. Note 
that h < h since the contact is deposited on top of the 
substrate. 



III. FORCES 

A. Electrostatic Force 

A numerical solution to the Poisson equation using a 
finite element method^ suggests that we can approxi- 
mate the capacitances of the junctions with a parallel 
plate capacitor model with an offset. The capacitances 
between the tube and the gate electrode C g and the tube 
and the drain electrode Cd are given bj^ 

Cd(z) = : x n C \^„ (6) 



C s (x) 



1- f(l -C /CV)' 



(7) 
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where Co is the drain junction capacitance at zero tube 
deflection, Ch is the drain junction capacitance when the 
tube is in mechanical contact with the drain electrode, 
and k = 0.5(3z^ — takes into account that the deflec- 
tion at the tube above the gate is smaller than at the tip. 
The value of Ch is estimated from experimental dataM 
and Co is calculated for parallel plate capacitors. The 
resulting capacitive forces are given by 




where qa, g are the charges on the drain and gate capaci- 
tors respectively, qd + Qg = 1, and q is assumed to be con- 
stant during an infinitesimal displacement of the tube. 

B. Van der Waals force 

The origin of the van der Waals (vdW) force is corre- 
lation of fluctuating higher order electrostatic moments 
between two charge distributions. The force is attractive 
and the interaction energy varies with separation as the 
inverse sixth power of separation in the range of separa- 
tions relevant to the nanorelay^ The total vdW-energy 
between the tube and the substrate plus contacts is calcu- 
lated in a continuum limit ^ The pairwise sum over the 
interaction between the individual atoms is transformed 
into a six-dimensional integral over the volumes of the 
interacting bodies. In this approximation we can write 
the vdW-force as 

-FvdW = — Vi?vdW 

X7 f f AVJV " g 6 (9) 

= - V / / dVldV 2 n 1 - rrrn 2 , 

Jv 1 Jv 2 |ri-r 2 | 6 

where V\ is the volume of the tube, V 2 is the volume of 
the contacts and substrate, and ni(ri) and n 2 (r2) are 
the atom densities in the two bodies. The values of the 
interaction parameter C§ depend on the species of the 
interacting atoms and is usually stated in terms of the 
Hamaker constant A — Tr 2 nin 2 Ce. The Hamaker con- 
stant for interaction between graphite-silicon dioxide is 
given by A — 0.07 aJ and for graphite-metal the corre- 
sponding value is 0.6 aJ»i£ 

Evaluating the integral in Eq. @ is time-consuming 
for an arbitrary geometry. To simplify this, we sep- 
arate -EVdW m two parts, one contribution from tube- 
electrodes interaction and one from tube-substrate inter- 
actions. The substrate is approximated by a semi-infinite 
plane and the drain electrode is assumed to be cylindri- 
cal and infinitely long. Furthermore, we assume that the 
interaction between the tube and the gate electrode is 
much smaller than the tube-drain interaction due to the 
much larger separation. Assuming the MWNT to con- 
sist of individual CNTs separated by 3.4 A, we transform 
the integral over one-dimension to sum over all individ- 
ual CNTs in the MWNT. Many of the integrals in © 



can be carried out analytically, leaving either one (for 
tube-substrate interaction) or three (for tube-electrode 
interaction) to be evaluated numerically. 

C. Short Range Forces 

The attractive part of the short range force F a tt origi- 
nates from a coupling between the tube and the contact 
when the electronic wave functions overlap. It can be 
roughly related to the tunneling conductance Gt = 1/ Rt 
between the tube and the contact asi& 

F m oc VG. (10) 

A more detailed analysis of the this force is rather com- 
plicated but does not change the qualitative features of 
the system's behavior^ When the core-core separation 
of the of the atoms in the tube and in the tip becomes 
small, however, the wave function overlap increases and 
we expect a strong repulsion due to the Pauli princi- 
ple. Taking this repulsive part F rop to be exponential 
we model the attractive and the repulsive short range 
forces by a Morse curved 

Fsr -Fatt ~t~ F rC p 

= ^ VeceTifcGoCe-^-^ - e"* <*-*>), ^ 

where / is a dimensionless shape factor of the order of 
one, ct and ec are the valence band widths of the tube 
and the drain electrodes, respectively, and Go = l/Ro is 
the conductance at contact. We define the location of the 
surface by F sr = 0. The core-core separation between the 
atoms in the tube and the drain electrode when the tube 
is in mechanical contact with the electrode is assumed 
to be 3 A. This offset is introduced to prevent the un- 
physical divergence in the vdW-energy at zero tube-drain 
electrode separation. 

IV. PHONON EXCITATION 

The only dissipative mechanism in the equation of mo- 
tion thus far is a viscous damping. Several other dissipa- 
tive mechanisms are present in a real system. One effect 
that is important for the device performance is the pos- 
sibility of phonon excitation in the drain contact when 
the tube bounces off the surface of the drain electrode. 
This surface dissipation is important for the switching 
behavior. 



A. Phonon model 

To estimate the dissipation associated with phonon 
emission we use a simplified classical one dimensional 
model of the atomic structure. The atomic lattice is 
modeled as a collection of atoms, each connected to the 
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Figure 2: The magnitude of the different forces in the system 
with parameters given in chapter V. The voltages V g = 7.5 V 
and V g = 0.01 V were used. Close to the contact located at 
x = h the forces are of the same order while electrostatic and 
elastic forces dominate far from the surface. 
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Figure 4: The remaining kinetic energy of the tube in units 
of incident kinetic energy as a function of incident velocity. 
Increasing the incident velocity increases the relative dissi- 
pation. The nanotube has an effective mass of 2.1 • 10 -21 
kg corresponding to the effective mass for the typical system 
described in chapter V. 
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Figure 3: Stability diagram with and without surface forces 
for the system with parameters given in Chapter V. The curve 
shows the positions of zero net force on the tube (or local equi- 
libria) as functions of gate voltage (at constant V s = 0.01 V) 
and deflection x (in units of h). The large arrows shows the 
direction of the force on each side of the curves, indicating 
one local equilibrium to be unstable. The required voltage 
for pulling the tube to the surface ( "pull-in voltage" ) is given 
by A (w 6.73 V). In agreement with reference Ha this volt- 
age is not significantly affected by surface forces. A tube 
at the surface will not leave the surface until the voltage is 
lower than the "release voltage", B and C in the figure. Note 
that A > B, C which indicate a hysteretic behavior in the 
IVg-characteristics, a feature significantly enhanced by sur- 
face forces. 



other through a harmonic potential, characterized by a 
spring constant K. Assuming linear dispersion relation 
we obtain the spring constan-fe2i 



where M is the mass of the atoms in the chain, a is the 
equilibrium separation between atoms, and c is the speed 
of sound in the electrode material. We consider the drain 
electrode to be made of gold, and thus, we use the lattice 
constant of gold a = 407.8 A, the mass of one gold atom 
M = 197 amu, and the sound velocity in gold c = 2000 
m/s. 

On the atomic scale the nanotube is a very large object 
and, thus, the tube will interact with the surface over an 
extended region while bouncing. We estimate that the 
impact area of the nanotube corresponds to about 50 
gold atoms (an impact area of the order of 2 nm 2 ) in the 
contact. Instead of considering the motion of each atom 
in the lattice, we consider the motion of layers of atoms, 
each layer made up of 50 atoms. The nanotube, charac- 
terized by the effective mass m e ff, impinges on the model 
lattice with an initial velocity v- ln , interacts with the top- 
most layer through a repulsive force F Tep and bounces. 
We use a total number of 100 layers, a number chosen 
such that the boundary conditions at the last layer have 
only minor importance. The motion of the lattice layers 
and the nanotube is governed by classical dynamics and 
the resulting one-dimensional problem is solved numeri- 
cally. 

The one-dimensional model we use does not, e.g., take 
into account phonons with different polarizations, hence, 
the model most likely underestimates surface dissipation. 
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V. RESULTS 
A. Effects of Surface Forces 

The surface forces introduce constraints on the design 
parameters. For the operation of the relay it is necessary 
that 

- kx + F c + F vdw + F sr > (13) 
-kx + F vdW + F sr < (14) 

in order to both enable pulling the tube to the contact 
(|13fl and for the tube to release when the voltage is turned 
off 1)1 4[l. In particular, the latter condition is hard to 
fulfill - a manifestation of the ubiquitous stiction problem 
in nano-science. 

We investigate the behavior of a typical system where 
Eqs. i|13|) and l|14fl are satisfied. The geometrical param- 
eters of this system are tube length L = 120 nm, inner 
tube diameter D\ = 2 nm, outer tube diameter D a = 8.8 
nm, tube contact separation at zero deflection h = 5 nm, 
relative position of the gate z r = 0.7, and Young's mod- 
ulus E = 1 TPa. With these geometrical parameters the 
effective mass is m c ff ~ 2.1 • 10 -21 kg. Other parameters 
are source junction impedance Z = 8 kf2, quality factor 
of the viscous damping Q = 250, capacitance at zero de- 
flection Co = 1.3 • 10 -19 F, and capacitance at contact 
C h = 5 ■ 10~ 18 F. 

The magnitudes of the forces in a typical system are 
depicted in Fig. [21 as a function of the deflection of the 
tube. When the tube tip is near the drain electrode, all 
forces are comparable, while far from contact the electro- 
static and mechanical forces dominate. 

The effects of the surface forces can be visualized by 
means of a stability diagram which shows the positions 
of zero net force (local equilibria) on the cantilever as a 
function of gate voltage and deflection. A typical stabil- 
ity curve for the system is given in Fig. [31 A stability 
diagram with more than one local equilibrium for a spe- 
cific voltage results in a hysteretic behavior in the IV g - 
characteristics (more than one position of zero net force 
means that there actually are three such positions, one of 
which is not stable) . This is because the net force is pos- 
itive to the right of the curve and negative to the left of 
the curve, so that, when lowering the voltage for a tube 
in contact with the drain electrode, it will not release un- 
til there is no stable position at the surface. Fig.[31shows 
that this effect is significantly increased by the surface 
forces. 

The "sticking problem" corresponds to a stable nan- 
otube position for zero voltage at the surface, i.e. fail- 
ure to satisfy l|14f) . This problem can be alleviated 
e.g. by choosing stiffer or shorter tubes. This in turn 
leads to higher voltages needed to satisfy equation (|13fl . 
With high voltages and small distances, the electric field 
strengths can be very high. Under these conditions field 
emission may be important. Simulations show, however, 
that design parameters can be chosen so that this has a 
negligible effect^ 
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Figure 5: Step response (gate voltage step (0 — > —7.5 V) 
applied at t — 2(j,s) of the system with parameters given in 
chapter V without phonon dissipation (elastic bounce) and 
with surface dissipation (inelastic bounce). The dissipative 
mechanism reduces the total time of switching by nearly two 
orders of magnitude. The noise reflects the stochastic nature 
of the tunneling between the tube and the drain electrode 
(shot noise). 
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Figure 6: The IV g -characteristics for a typical system with 
the parameters given in Chapter V has a large hysteresis loop. 
The arrows indicate up- or down sweep of V s (at constant 
V s = 0.01 V). The system switches to the conducting state at 
V g ~ 6.7 V which agrees well with the expected value deduced 
from Fig. 01 The reverse transition when lowering the gate 
voltage takes place at V g « 3.1 V which is slightly larger than 
the expected value. This small discrepancy is due to small 
vibrations of the tube and non-adiabatic voltage sweep. 
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(a) T w = 0.015 ns. The 
system stays in the 
conducting 1-statc. 



(b) T w = 0.020 ns. The 
system switches to the 
non-conducting 0-state. 



(c) T w = 0.7 ns. The (d) T w = 0.8 ns. The 

system stays in the system switches the 

non-conducting 0-state. conducting 1-state. 

Figure 7: The system's response for gate voltage pulses of 
duration T w applied at t = 0.2 /is when biased within the 
hysteresis loop (V B = 5 V). The parameters of the system are 
given in chapter V. (a,b) low voltage pulse V s (t) — V to 
write a '0'; (c,d) high voltage pulse V g (t) — 7.5 V to write a 
'T. Note the asymmetry between the write times of a '0' and 
aT. 



B. Effects of phonon dissipation 

As the incident velocity of the tube increases a larger 
percentage of energy is transferred from the tube to the 
substrate. For our particular set of simulation parame- 
ters we get a loss function according to Fig. For high 
incident velocities we get high dissipation whereas for low 
velocities the relative dissipation gets smaller. This prop- 
erty makes relaxation towards the surface fast even for 
high initial velocities. 

The main effect of exciting phonons in the drain con- 
tact is to change the switching dynamics of the relay. 
Without this dissipative mechanism, relaxation of the 
tube towards a stationary current carrying state is slow 
due to the fact that the tube bounces off the contact 
surface many times before coming to a rest. This time 
scale is significantly decreased when phonon dissipation is 
included in the model. A comparison between a simula- 
tion showing the response to an applied step gate voltage 
without dissipation and an equivalent simulation with 
surface dissipation, depicted in Fig. [5] show that the time 
scale of switching is decreased by two orders of magni- 



tude. The reverse transition is not affected by including 
phonon dissipation in the model. 



C. Memory Element Application 

We have performed numerical simulations of the sys- 
tem described above and the resulting I V g -characteristics 
at a constant V s — 0.01 V are shown in Pig. El The hys- 
teresis is large, which makes a memory element a promis- 
ing application of the device. In such a memory element 
the conducting state can be defined as a logical 1 and 
the non-conducting state as a logical 0. The system is 
biased at a gate voltage chosen to lie within the hys- 
teresis loop where the system is bistable. The element 
is written by applying a high/low voltage for a specific 
write time, T w . The low voltage is typically zero and high 
voltage is chosen such that the condition (|13fl is fulfilled. 
The write times must be long enough to ensure that the 
correct state is reached, and depend on the chosen write 
voltages. 

We have estimated the write times of the memory ele- 
ment by numerical simulations. Such simulations for the 
chosen set of parameters are depicted in Fig. [7] The write 
time for — > 1 is approximately 0.8 ns and for 1 — * 
the write time is about 0.02 ns. The difference between 
these time scales can be explained by comparing the two 
write processes: When the tube moves from a stationary 
1, the tunneling resistance increases exponentially with 
tube contact separation and the current will stop soon 
after the voltage is switched off. In the — > 1 transition 
the tube bends relatively quickly to the drain electrode, 
but tends to bounce off the surface. The kinetic energy 
of the mechanical motion must be dissipated before the 
tube relaxes toward the stationary 1-state. This relax- 
ation is enhanced due to phonon dissipation in the drain, 
but we cannot expect it to remove the asymmetry com- 
pletely. The time scale corresponding to the — > 1 tran- 
sition sets the limit for maximum operating frequency, 
approximately 1 GHz for the parameters we have used. 



VI. CONCLUSIONS 

We have incorporated short range and vdW-forces into 
a model of a three terminal nanorelay. The main ef- 
fect of these forces is to increase the hysteresis in the 
IV^-characteristics making a memory element a promis- 
ing application of the relay. We have investigated the 
switching characteristics for such a memory element and 
conclude that the write dynamics is limited by the — > 1 
transition. The transition dynamics of the — ► 1 tran- 
sition depends to a great extent on the dissipative sur- 
face processes when the tube bounces, which reduce the 
transition time by two orders of magnitude, whereas the 
1 — > transition is unaffected by such processes. 
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